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Abstract 

An algorithm is presented allowing the construction of fast Fourier 
transforms for any solvable group on a classical computer. The 
special structure of the recursion formula being the core of this 
algorithm makes it a good starting point to obtain systematically 
fast Fourier transforms for solvable groups on a quantum computer. 
The inherent structure of the Hilbert space imposed by the qubit 
architecture suggests to consider groups of order 2 n first (where n is 
the number of qubits). As an example, fast quantum Fourier trans- 
forms for all 4 classes of non-abelian 2-groups with cyclic normal 
subgroup of index 2 are explicitly constructed in terms of quantum 
circuits. The (quantum) complexity of the Fourier transform for 
these groups of size 2 n is 0(n 2 ) in all cases. 
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1 Introduction 



Quantum algorithms are a recent subject and of possibly central importance in 
physics and computer science. It has been shown that there are problems on 
which a putative quantum computer could outperform every classical computer. 
A striking example is Shor's factoring algorithm (see p2| ). 

Here we adress a problem used as a subroutine in almost all known quantum 
algorithms: The quantum Fourier transform (QFT) and its generalization to 
arbitrary finite groups. 

In classical computation there exist elaborate methods for the construction of 
Fourier transforms (see Beth ||, Q, Clausen ||, ||, and Diaconis/Rockmore [[|) 
so it is highly interesting to adapt and modify these methods to get a quantum 
algorithm with a much better performance (with respect to the common quantum 
complexity model) as on a classical computer. 

First attempts in this direction have been proposed by Beals [0 and H0yer 



T2| . In this paper we present an algebraic approach using representation theory 
which can be seen as a first step towards the realization of a large class of QFTs 
on a quantum computer. 



2 Generalized Fourier Transforms 

Fourier transforms for finite groups are an interesting and well studied topic 
for classical computers. We refer to 0, ||, |15| , fT9|| as representatives for a 



vast number of publications. The reader not familiar with the standard notation 
concerning group representations should refer to these publications or to standard 
references as |§ or lf2~lH . 

For the convenience of the reader we first recall the definition of the generalized 
Fourier transforms for a finite group G and explain the representation theoretical 
point of view we are going to take. 

Each isomorphism 

k 

$ : CG »0C rf ' xdl 



i=l 



between the group algebra of G and the Wedderburn components is called a 
Fourier transform for the group G. A particular isomorphism is fixed by picking 
a system pi, . . . , pk of representatives of irreducible representations of G and 
defining $ as the linear extension of the mapping g i— > 0f =1 Pi(g), g £ G (of 
course deg(pj) = di). Any "fast" algorithm for the evaluation of $ is called a fast 
Fourier transform for G. 

In order to obtain a matrix representation for the linear mapping $ one usu- 
ally fixes natural bases L in CG and V in 0f =1 C diXdi . This is done by choosing 
an ordering L = (gi, . . . ,g\c\) on the elements of G and an ordering V on the 
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elementary matrices E k j (1 at position (k, I), else) which correspond to the coef- 
ficients of the irreducible representations appearing in the Wedderburn decompo- 
sition. The Fourier transform $ then is represented by a matrix Ml^i e C' G ' X ' G '. 
Since Ml,l' is a base change between two orthonormal bases (with respect to the 
standard hermitean scalar product on CG) it must be unitary. 

The matrix M LjL > can also be characterized by another property. Let be the 
regular representation obtained by right multiplication of G on L. Conjugating 
with Ml^i yields (up to a permutation matrix P) a direct sum of irreducible 
representations with the property that equivalent irreducibles are equal, i.e. 



On the other hand every matrix with this property corresponds to a Fourier 
transform with respect to natural bases. 

As an example let G = Z n = (x \ x n = 1) be the cyclic group of order n with ir- 
reducible representations pi = (x i— > u l n ), i — . . .n — 1, where u n denotes a prim- 
itive n-th root of unity. With respect to the natural bases L — (x l \ i — . . .n — 1) 



and L< = (E 1A , E n>n ) the matrix M LtI/ = \ \ i, j = . . .n - 1] = DFT 



is the discrete Fourier transform well-known from signal processing. 

We will refer to the notion of a fast Fourier transform as a fast algorithm 
for the multiplication with M LjL i. Of course, the term "fast" depends on the 
chosen complexity model. Since we are primarily interested in the realization of 
a fast Fourier transform on a quantum computer (QFT) we first have to define 
the measure of complexity on this architecture. 



In this paper we think of a quantum computer to consist of a quantum register 
which in turn consists of n qubits each of which provides a 2-dimensional Hilbert 
space. Thus the possible operations this computer can perform are given by the 
unitary group U(2 n ). 

To study the complexity of unitary operators on n-qubit quantum systems we 
introduce the following two types of building blocks: 

• Local unitary operations on the qubit % are of the form 

C/W ;= 1 2 ;_! (g) U® l 2 n-i, 

where U is an element of the unitary group 14(2) of 2 x 2-matrices. 

• The controlled NOT gate (also called measurement gate) between the qubits 
i (control) and j (target) is defined by 



<f) M L,L> = ( Pl e . . . e p k ) p fulfilling Pi P j =>■ Pi = Pj . 



3 The Complexity Model 
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CNOT^ 



Figure 1: Elementary quantum gates 

when restricted to the tensor component of the Hilbert space which is 
spanned by the qubits i and j. Note that the controlled not is nothing 
but an XOR on the resp. components. 

In the graphical notation using quantum wires these transforms are written as 
shown in figure The lines represent the qubits and the least significant bit is 
the lowest. We assume that these so-called elementary quantum gates can be 
performed with cost 0(1). 

These two types of gates suffice to generate all unitary transformations, which 
is the content of the following theorem from [|J. 

Theorem 3.1 The set Q = {17«, CNOT^ | U G U(2), i,j = l...n, i ^ j} is 

a generating set for the unitary group U(2 n ). 

This means that for each U G U{2 n ) there is a word vj\Wi ■ ■ ■ Wk (where Wi G Q 
for % = 1 . . . k is an elementary gate) such that U factorizes as U — W1W2 ■ ■ ■ Wk- 

In general only exponential upper bounds for the minimal length occuring in 
factorizations have been proved (see |1[) but there are many interesting classes 
of unitary matrices in U(2 n ) affording only polylogarithmic word length, which 
means, that the minimal length k is asymptotically 0(p(n)) where p is a polyno- 
mial. 

In the following we give examples of some particular unitary transforms ad- 
mitting short factorizations which will be useful in the rest of the paper. 

• The symmetric group S n is embedded in U{2 n ) by the natural operation of 
S n on the tensor components (qubits). Let r G S n and Il r the correspond- 
ing permutation matrix on 2 n points. Then I1 T has a 0(ri) factorization 
as shown in [TS|]. As an example in figure |2] the permutation (1,3,2) of 
the qubits (which corresponds to the permutation (2, 5, 3) (4, 6, 7) on the 
register) is factored according to (1,3,2) = (1,2)(2,3). 

• Following the notation of [|l|] we denote a fc-times controlled U by Ak(U). 
Lemma 7.2 and lemma 7.5 in [Q] show that for U G U{2) the gate A n _i(/7) 
can be realized with gate complexity 0(n 2 ) and Ak(U) with 0(n) for k < 
n — 1 on n qubits. 
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Figure 2: Factorization (1,3,2) = (1,2)(2,3) 

The Fourier transform DFT 2 " can be performed in 0(n 2 ) elementary oper- 
ations on a quantum computer (see p2fl , J7|). 



Let P n G be the cyclic shift permutation which acts on the states of 
the quantum register as x i— > x + 1 mod 2™. Obviously the corresponding 
permutation matrix is the 2 n -cycle 



/0 1 
1 



V 1 



\ 



o 1 



A factorization of P n as a product of multiple controlled NOTs as shown in 
figure |3] needs 0{n 2 ) basic operations^. 
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Figure 3: Quantum circuit for the 2 n -cycle P n 



• Let U G U{2 n ). The cost for a (single) controlled [/ is settled by the 
following lemma. 



Lemma 3.2 IfU G U(2 n ) can be realized inO(p(ri)) elementary operations 
then Ai(U) G U(2 n+1 ) can also be realized in 0(p(n)) basic operations. 

lr This quantum circuit, similar to the classical carry-look-ahead logic, has been found by 
Markus Grassl following discussions with Amir Fijany in May 1998 
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Proof: First we assume without loss of generality that U is written in 
elementary gates. Therefore we have to show that a doubly controlled 
NOT and a single controlled U G W(2) can be realized with a constant 
increase of length. This follows from |I| . □ 



4 Creating Fast Fourier Transforms 

In section |2| we have explained that calculating a Fourier transform for a group G 
is the same as decomposing a regular representation of G with a matrix A into 
irreducibles pi up to a permutation P such that equivalent irreducible summands 
are equal, i.e. 

(f) A = A' 1 ■ <p- A= ( Pl ® ...® Pk ) P fulfilling Pi = Pi = Pj . 

A "fast" Fourier transform (on a classical computer) is given by a factoriza- 
tion of A into a product of sparse (w.r.t. the architectural complexity measure) 
matrices. 

In this section we present (without proof) a number of theorems and lemmata 
yielding an algorithm to calculate fast Fourier transforms for any solvable group G 
(on a classical computer). The same algorithm serves as a good starting point to 
obtain quantum Fourier transforms if we assume the "quantum wires" to possess 
a suitable number of states. 



The statements in this section all are taken from the first chapter of fl9 
where decomposition matrices and constructive representation theory in general 
is investigated. There the objective is the construction of decomposition matrices 
for monomial representations which can be viewed as a generalization of Fourier 
transforms. 

The following theorem provides the crucial formula needed to obtain a fast 
Fourier transform of G by decomposing a regular representation stepwise along 
a composition series of G. The formula has been known (see ||) to yield fast 
Fourier transforms for solvable groups on a classical computer (counting addi- 
tions and multiplications). Their tensor structure, however, also fits well to the 
special architecture of a quantum computer. The general constructive form of 



the following theorem as presented is due to [[HJ where furthermore an improved 
recursion formula for the classical architecture can be found. 

We use the following convention for the induction of a representation of 
H < G with transversal (i.e. a system of representatives for the right cosets) 
T — (ti, . . . , t n ): 



Tr G){x) = (f){t i xt j x ) | i, j 



. n 



where <j)(y) = <j)(y) for y G H and the all-zero matrix else. A regular represen- 
tation cf) is given by an induction <p = (lg j> G) where 1e denotes the trivial 
representation of the trivial subgroup E < G. 
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Theorem 4.1 Let N <j G a normal subgroup of prime index p with (cyclic) 
transversal T = (t°, t 1 , . . . , t^ p ~^) and <p a representation of degree d of N which 
has an extension <p to G (see figure Suppose that A is matrix decomposing <f> 
into irreducibles, i.e. (fi A = p = p 1 © . . . © p k and that p is an extension of p to 
G. Then 

p-i 

B = (l p ®A)-D-(DFT p ®l d ), where D = ($)p(t)\ 

1=0 

is a decomposition matrix for G, more precisely 

(0 Tt G) B = ©A r p, 

where \ : t t— > bf , % — 0. . .p — 1, are the p 1-dimensional representations of G 
arising from the factor group G/N . 



G 



ext 



ext 



N P 

dec 



Figure 4: Situation in theorem W7L 



In the case of an abelian group G the formula yields exactly the well-known 
Cooley-Tukey decomposition, hence D can be viewed as a generalized Twiddle 
factor. In the case of G being a direct product G = N x G/N the Twiddle matrix 
D vanishes. Since we want to apply theorem [TT] to a regular representation we 
need the following lemma. 

Lemma 4.2 Let N <j G a normal subgroup of prime index p and 4> any regular 
representation of N. Then <p (and hence all of its conjugates) has an extension 
(f) to G. Furthermore = 0* for all t G G. (<$" : ih cj){txt~ l ) is called the inner 
conjugate of <fi by t). 



In order to obtain an algorithm from theorem |4.1| we are faced with two problems. 
The first is the calculation of the Twiddle matrix D which is essentially the 
problem of extending p to p and evaluating it at t. Suppose we are given A 
and p which is a direct sum of irreducibles, p — p\ © . . . ffi p^, with equivalent 
summands being equal. Because of lemma |4.2| G/N operates on the irreducibles 
Pi via inner conjugation (explained in the lemma above). According to Clifford's 
Theorem (see e.g. B, pp. 88) exactly one of the following two cases applies 
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to each summand p«: Either pi = p\ and pi can be extended to G or pi ^ p\ 
and pi ]t G is irreducible. In the first case the extension may be calculated by 
Minkwitz' formula (see [fT/H), in the latter case the direct sum pi® p\@ ■ ■ - ® pf P ^ 
can be extended to G by pi ]t G. We do not state Minkwitz' formula here since 
we will not need it in the special cases treated in section ||. 

The second problem arises from the fact that the decomposition Xi ■ ~p in 



theorem |4.1| does not satisfy the property of equivalent summands being equal. 



This can be achieved using the following lemma. 

Lemma 4.3 Let N <j G a normal subgroup of prime index p with transversal 
T = t 1 , . . . ,*Cp-i)). Suppose that p is an irreducible representation of degree 
d of N satisfying p ^ p l and Xi : t i— > uj % p is an irreducible representation of G 
arising from G / N . Then 

(\ ■ (P Tt G)) mid = P 1tG, D = diag(l, u p ,..., uWy. 

Now we are ready to formulate the algorithm which constructs a fast Fourier 
transform for G from a fast Fourier transform of a normal subgroup of prime 
index. 

Algorithm 4.4 Let N < G a normal subgroup of prime index p with transversal 
T = (t°,tV--,* (p ~ 1) )- Suppose that is a regular representation of N with 
decomposition matrix A: 



lA 



pi © • • • © Pk fulfilling pi = pj =>- pi = pj. 



A decomposition matrix B for the regular representation <f> | T G can be obtained 
as follows. 

1. Determine a permutation matrix P rearranging the pi, i = l...k, such 
that the extendable pi (i.e. those satisfying pi = p\) come first followed by 
the others ordered into sequences of length p equivalent to pi, p\, . . . , pf P 11 . 
(Note: These sequences we need to equal p^ p\, . . . , pf P 1} which is estab- 
lished in the next step). 

2. Calculate a matrix M which is the identity on the extendables and conju- 
gates the sequences of length p to make them equal to pi, pj, . . . , pf P x) . 

3. Note that A-P-M is a decomposition matrix for cf>, too, and let p = <p A ' p ' M . 
Extend p to G summandwise. For the extendable summands use Minkwitz' 
formula, the sequences pi, p\, . . . , pf P x) can be extended by pi |t G as stated 
above. 

p-1 

4- Evaluate p at t and build D = @ p(t) 1 - 

i=0 
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Figure 5: Coarse Quantum circuit visualizing algorithm ^L4 



5. Construct a blockdiagonal matrix C with lemma conjugating 0fr o Aj • p 
such that equivalent irreducibles are equal. C is the identity on the extended 
summands. 



Then 



B = (l p ® A ■ P ■ M) ■ D ■ (DFT 

is a decomposition matrix for <p f> G. 



C 



(1) 
□ 



It is obviously possible to construct fast Fourier transforms on a classical com- 
puter for any solvable group by recursive use of this algorithm. Note that the 
consideration of T- adapted representations (see 0) here is unnecessary: The 
irreducibles are constructed along with the decomposition matrices. 

Since we restrict ourselves to the case of a quantum computer consisting of 
qubits, i.e. two-level systems, we apply algorithm |4.4j to obtain QFTs for 2-groups 
(size is a 2-power). In this case the two tensor products occuring in (JJ) fit very 
well to yield a coarse factorization as shown in figure The remaining problem, 
however, is the realization of the matrices A, P, M, D, C in terms of elementary 
building blocks as presented in section |3|. At present this realization remains 
a creative process which might be performed by hand if an arbitrary class of 
groups is given. In section |5] we will apply algorithm iA to a class of non-abelian 
2-groups. 



5 Generating QFTs for a class of 2-groups 

In the case of G being an abelian 2-group the realization of a fast quantum 
Fourier transform has been settled by |T3 |. Clearly, this case is covered by the 



method presented here (see the remarks following theorem 4~1). In this section we 
will apply algorithm [4.4| to the class of non-abelian 2-groups containing a cyclic 
normal subgroup of index 2. Fast quantum Fourier transforms for these groups 
have already been constructed by H0yer in [O. 
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According to HT3" |, p. 90/91 there are for n > 3 exactly four isomorphism types 



of non-abelian groups of order 2 n+1 affording a cyclic normal subgroup of order 
2 n : 

(i) The dihedral group D 2 n+i — (x, y | x 2 ™ = y 2 — 1, a; 2 ' = x -1 ). 

(ii) The quaternion group Q 2 n+i — (x, y \ x 2 " = y 4 = 1, x y = x^ 1 ). 
(hi) The group QP 2 n+i = (x,y \ x 2 " = y 2 = 1, x y = x 2 " 

(iv) The quasidihedral group QD 2n +i = (x, y \ x 2 " = y 2 = 1, x y = x 2 " 

Observe that the extensions (i), (iii), and (iv) of the cyclic subgroup Z 2n = (x) 
split, i.e. the groups have the structure of a semidirect product of Z 2 n by Z 2 . 
The three isomorphism types correspond to the three different embeddings of 
Z 2 = (y) into (Z 2 n) x ^ Z 2 x Z 2 n- 2 . 

5.1 QFT for the dihedral groups D 2 n+i 

In this section we construct a QFT for the dihedral groups D 2 n+i step by step 



according to algorithm ^] and explicitly state the occuring quantum circuits. 

Let G = D 2 n+i = (x,y \ x 2 ™ = y 2 = 1, x y = x~ l ) with normal subgroup 
N = (x) < G of index 2 and transversal T = (l,y). We consider the regular 
representation = (1# AT) | r G of G with S* = (l,x, . . . ,x 2 " -1 ). Obviously 
the regular representation (1# js AT) of A" is decomposed by A = DFT 2 « into 
Po © • • • © P2"-i where pi = (x h- > cu^O- Now we are ready to apply algorithm [14] 
to obtain a decomposition matrix S for 0. For convenience we denote u 2 n simply 
as uj and the Hadamard matrix as 

1. Since pf (x) = piiyxy' 1 ) = p^a; -1 ) = p 2 n_i(x) we see that there are exactly 
two extendable pi namely for i = 0, 2 n ~ 1 . The sequences of inner conjugates 
are given by pi, p 2 ^-i, i ^ 0, 2 n_1 . We need a permutation P reordering the 

Pi as 

PO, P2™- 1 ) Pi) P2™-1, • • • , Pi, P2"-i, • • • , P2™- 1 -!? P2"- 1 +l • 
extendables pairs of inner conjugates 

This can be accomplished by the circuit given in hgure |^ since the n-cycle 
on the qubits which is performed first yields a decimation by two on the 
indices, i.e. the indices 0, . . . ,2 n_1 — 1 have found their correct position. 
The only thing which remains to do is to perform the operation x 1— > — x 
on the odd positions. This can be done by an inversion of all (odd) bits 
followed by a x 1— > x + 1 shift P n -i on the odd states of the register. 
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Figure 6: Ordering the irreducibles of Z 2 n <j D 2 n+i 

2. M can be omitted since all the are of degree 1. 

3. Let 4> A ' P = p. We extend p summandwise to p: 

• po = In can be extended by 1q. 

• p 2 n-i can be extended through p 2 n-i(y) = 1. 

• The sequences pi © p2«-i, i 7^ 0, 2 n ~ 1 can be extended by pi f t G: 



Pi |r G : x I u -i J, y 



1 

1 



4. Evaluation of p at the transversal T yields the Twiddle matrix 



D = p(l)©p(y) 
/I , 



1 



V 



I 1 / 



D is realized by the quantum circuit given in figure [7|. 
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Figure 7: Twiddle matrix for D 2 n+i 
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Figure 8: Equalizing inductions 
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Figure 9: Complete QFT circuit for the dihedral group D 2 n-\ 



5. According to lemma (O] the matrix C has the following diagonal form: 

C= l 2 n ©diag(l, 1,1,-1,..., 1,-1), 

V v ' 

2«-l-l pairs 

which is realized by the quantum circuit given in figure ||. 

Summarizing we obtain that 

B = (l p ® A ■ P ■ M) ■ D ■ (DFT P <g> 1,^) • C 

is a decomposition matrix for and a fast quantum Fourier transform for G. The 
whole circuit is shown in figure || 



5.2 QFT for the groups Q 2 n+i, QP2"+i, and QD 2 n+i 

In the following we give the circuits for the groups Q%n+i, QP2n+i, and QD 2 «+i. 
In all cases we have (x) = N <j G so that algorithm [4.4| has to be performed only 
once for the last step. For the sake of brevity we will state only those parts of 
the circuit which differ from the dihedral group. Some of the essential properties 
of the groups are summarized in tabular [TJ. We use the same notation as in the 
last section. 



Group 


Inner conjugates of Z 2 n 


No. of 1-dim 
irreducibles 


No. of 2-dim 
irreducibles 


D2 n + 1 


Pi, P2 n -i 


4 


2 n-l _ 1 




Pi, P2 n -i 


4 


2 n-l _ l 




Pi, Pj(2 Il - 1 +l)mod2" 


2 n 


2 n-2 


QD211+1 


Pi, Pi(2™- 1 -l)mod2" 


4 


2 n-l _ X 



Table 1: A class of non-abelian 2-groups 
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Figure 10: Twiddle matrix for Q 2 n+i 



Figure 11: Permutation for QP 2 n+i 



Q 2 n+i: The irreducibles pi extend or induce in the same way as in the 
dihedral case. Hence the QFT only differs in the Twiddle matrix D since 
for a not extendable pi we have 



(pi ] T G){y)=\ 

Thus the Twiddle matrix D is given by 

/I , 



1 
-1 



D 



-1 



V 



-I 1 / 



and can be realized by the circuit in figure |T0 

detern 
Hence 



QP 2 n+i: To determine which pi are extendable we use p\{x) = Pi(yxy x ) 

t ,2"- 1 +l^ 



pi = p\ u =u 



j.(2"-l+l) 



& U;^- 1 = 1 & 2 I i 



and there are exactly 2™ 1 extendable pj. The reordering permutation P 
has the easy form shown in figure |TT], and the matrix D is given by 



D 



low- 



/l 1 

V 



I 1 / 



which is simply a doubly controlled not as visualized in figure 
The matrix C then is given by figure 

QD 2 n+i: Here we have pf(x) = pi(x 2 ™ and 



p . = p y^^ = ^•(2- 1 -i) ^ w i-(2»--2) = 1 ^ • = a 2 



n-1 
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Figure 12: 
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Figure 13: Equalizing for QP 2 
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Figure 14: The permutation for the QD 2 ^ 



Thus everything is the same as in the dihedral case beside the orderin 
permutation P which takes the more complicate form shown in figure [14. 



Investigation of the quantum circuits yields the following theorem. 

Theorem 5.1 The Fourier transforms for the groups G — D 2 n^ Qi n i QP2 n > 
and QD 2 n can be performed on a quantum computer in 0(log 2 |G|) elementary 
operations. 



Proof: We can treat the four series uniformly, since the Fourier transforms all 
have the same decomposition pattern. First, in all cases a Fourier transform for 
the normal subgroup Z 2 n-i is performed with cost of 0(n 2 ) basic operations. The 
reordering permutation P, the Twiddle matrix D, and the equalizing matrix C 
cost 0{n 2 ) in case of D 2n , Q 2 ™, and QD 2 n due to lemma and example [| For 
QP 2 n we need only 0(1) operations for P, D, and C . □ 

All presented Fourier transforms have been implemented by the authors in the 
language GAP |2(J using the package AREP [|ll[] which will be available soon as 
a GAP share package. 
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6 Conclusions and Outlook 



A constructive algorithm has been presented allowing to attack the problem of 
constructing fast Fourier transforms for 2-groups G on a quantum computer built 
up from qubits. For a certain class of non-abelian 2-groups the algorithm has 
been successfully applied. All the QFTs created are of computational complexity 
0(log 2 \G\) like in the case of the cyclic group Z 2 ™. The main problem imposed 
by the implementation of certain permutation and block diagonal matrices has 
been solved efficiently 

Using the recursion formula from theorem [O it should be possible to con- 



struct QFTs for other classes of groups as well as to realize certain signal trans- 
forms on a quantum computer by means of symmetry-based decomposition (see 
0,0,0). 



We are indebted to Markus Grassl for helpful comments and discussions. Part 
of this work was presented and completed during the 1998 Elsag-Bailey - I.S.I. 
Foundation research meeting on quantum computation. 
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